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Abstract 

Previous measurements of i? = a{e^e~ 
hadrons)/(T(e"'"e~ — > at high are reviewed. 

Recent R measurement results, including those from the 
Beijing Spectrometer Experiment, are described. The 
present status of R measurements and future measurement 
possibiUties are summarized. 

1 INTRODUCTION 



Here, R measurements above the center-of-mass (cm) 
energy corresponding to the J/il) are reviewed. R mea- 
surements at lower energy are described in ref. ^ ^. Since 
the theoretical = ((; — 2)/2 is more dependent on low 
energy than on high energy R values, we will concern our- 
selves here primarily with the effect of R measurements 
on a{q^). In the following sections, previous R measure- 
ments will be summarized, recent measurements by the 
Beijing Spectrometer Experiment (BES) and CLEO will be 
described, and the current R status will be reviewed. 
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The QED running coupling constant evaluated at the Z 
pole, a(M^), and the anomalous magnetic moment of the 
muon, Ufj, = {g — 2)/2, are two fundamental quantities 
that are used to test the Standard Model (SM). The dom- 
inant uncertainties in both Q!(A/|) and a^^' are due to 
the effects of hadronic vacuum polarization, which can- 
not be reliably calculated in the low energy region. In- 
stead, with the application of dispersion relations, experi- 
mentally measured R values are used to determine the vac- 
uum polarization, where R is the lowest order cross sec- 
tion for e+e^ ^ 7* ^ hadrons in units of the lowest- 
order QED cross section for e+e^ l^-^f^^^ namely 
R = cr(e+e^ hadrons)/cr(e+e^ where 
cr(e+e~ = 0-°^ = 47ra^(0)/3s. Improved pre- 

cision for a(Mf ) also narrows the allowed range of the 
Higgs mass prediction using electro weak loop corrections. 
For much more detail on the importance of R measure- 
ments, see references |^, |^]. 

The uncertainty in a(M^) is introduced when it is ex- 
trapolated to the Z-pole. 
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where ao is the fine structure constant, which is known 
very precisely, and Aa{q^) is the vacuum polarization 
term. 

Aa{q') - Aaiiq') + K'IM<i') + Atopa(<z') 

The leptonic vacuum polarization, Aa/(g^), can be calcu- 
lated theoretically, and the top contribution, AtopQ!(g'^) is 
absorbed as a parameter in the SM fit. The dominant un- 
certainty is then due to the effects of hadronic vacuum po- 
larization, Aj^'^^, which is calculated with experimentally 
determined R values using dispersion relations. 



2 PREVIOUS R MEASUREMENTS 

Figure |l] shows the summary of R measurements as sum- 
marized by the Particle Data Group [|[]. These are selected 
measurements; older measurements with even larger errors 
are not shown. Also systematic normalization errors (5 - 
20 %) are not included on the points shown. 
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Figure 1 : Summary of R values from PDG2000 [||] . Some 
older measurements with even larger errors are not shown. 
Systematic normalization errors (5 - 20 %) are not included 
on the points shown. 

We will be primarily interested in the center of mass 
(cm) energy region nij/^ < i?cm < 12 GeV. For higher 
energies, perturbative QCD (PQCD) can be used to de- 
scribe the behavior of i? as a function of energy . Anal- 
yses that use data below ~ 12 GeV are "data driven" ap- 



proaches "theory driven" approaches use PQCD to go 
much lower in energy [^. Ahhough less dependent on the 
quality of experimental data, the latter must make addi- 
tional theoretical assumptions. 

Figure || shows a 1995 summary plot of the cm energy 
region below 10 GeV [^. The R values used as a function 
of energy to evaluate the dispersion integral are indicated 
by the smooth line through the points, and the R value un- 
certainties ascribed in the various energy regions are indi- 
cated by the bands. In the 1 to 5 GeV region, the R uncer- 
tainty is taken to be 15 %. Note that the Crystal Ball data 
used in the 5 to 7 GeV region is unpublished. 
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Figure 2: Summary of 1995 R values from ref. [^. The 
uncertainties in the various energy regions are shown by 
the bands and given at the bottom of the figure. 



3 BES R MEASUREMENT 

Recently a detailed R scan was completed by the upgraded 
Beijing Spectrometer (BES-II) Experiment [p^. The anal- 
ysis is described in some detail here to demonstrate the 
complexity of R measurements. BESII is a general purpose 
solenoidal detector located at the Beijing Electron Positron 
Collider (BEPC), which is the only facility operating in the 
CM energy range from 2 to 5 GeV. The luminosity at the 
J/^ is ^ 5 X 10'^° cm^^s^^. BESII is described in detail 
elsewhere [pTjl. 

Experimentally, the value of R is determined from the 
number of observed hadronic events, -/V^^^, by the relation 
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where N^g is the number of beam-associated background 
events; J2i-^u^ = s,fi,T) are the numbers of lepton- 
pair events from one-photon processes and N^^ the num- 
ber of two-photon process events that are misidentified as 
hadronic events; L is the integrated luminosity; 6 is the 
radiative correction; Chad is the detection efficiency for 
hadronic events; and etrg is the trigger efficiency. 

In 1998 BES performed an initial measurement of R us- 
ing six scan points: 2.6, 3.2, 3.4, 3.55, 4.6, and 5.0 GeV 
[ |l2| ] . At each point separated beam runs were done for the 
study of beam gas background. In 1999, BES measured 85 
scan points with 1000 hadronic events per point. Sep- 
arated beam running was done at 24 energy points, while 
single beam running was done at 7 points [|lO|. 

3. 1 Event Selection 

The main sources of background for this measurement are 
cosmic rays, lepton pair production, two-photon processes 
and single-beam associated processes. Clear Bhabha 
events are first rejected. Then the hadronic events are se- 
lected based on charged track information. Special atten- 
tion is paid to two-prong events, where cosmic ray and lep- 
ton pair backgrounds are especially severe, and additional 
requirements are imposed to provide extra background re- 
jection [p^. 

An acceptable charged track must be in the polar angle 
region | cos9\ < 0.84, have a good helix fit, and not be 
clearly identified as an electron or muon. The distance of 
closest approach to the beam axis must be less than 2 cm 
in the transverse plane, and must occur at a point along the 
beam axis for which \z\ < 18 cm. In addition, the follow- 
ing criteria must be satisfied: (i) p < pbeam + 5 x cTp, where 
p and Pbeam are the track and incident beam momenta, 
respectively, and ap is the momentum uncertainty for a 
charged track for which p Pbear,i\ (ii) E < 0.6Ei,eani, 
where E is the barrel shower counter (BSC) energy asso- 
ciated with the track, and Ebeam is the beam energy; (iii) 
2 < i < ij, + 5 X CTj (in ns.), where t is the measured 
time-of-flight for the track, tp is the time-of-flight calcu- 
lated assigning the proton mass to the track, and at is the 
resolution of the barrel time-of-flight system. 

After track selection, event selection requires the pres- 
ence of at least two charged tracks, of which at least one 
satisfies all of the criteria listed above. In addition, the to- 
tal energy deposited in the BSC (Esum) must be greater 
than 0.28Eiieam, and the selected tracks must not all point 
into the forward (cos 6* > 0) or the backward (cos 9 < 0) 
hemisphere. 

For two-prong events, residual cosmic ray and lepton 
pair (e+e^ and backgrounds are removed by re- 

quiring that the tracks not be back-to-back, and that there 
be at least two isolated energy clusters in the BSC with 
E > 100 MeV that are at least 15° in azimuth from the 
closest charged track. This last requirement rejects radia- 
tive Bhabha events. 

These requirements eliminate virtually all cosmic rays 



and most of the lepton pair (e+e~ and events. The 

remaining background contributions due to lepton pairs 
(Nil), including t+t^ production above threshold, and 
two-photon events (N^-y) are estimated using Monte Carlo 
simulations and subtracted as indicated in Eq. (1). 

Cuts used for selecting hadronic events were varied in 
a wide range, e.g. cos 9 from 0.75 to 0.90, Egum from 
Q.lAEbeam to Q.32Ei,eam, to estimate the systematic error 
arising from the event selection, which turned out to be the 
dominant systematic error as indicated in Table |[ 

3. 2 Beam Associated Background 

The largest background is due to beam associated back- 
ground, Nbg. To determine the level of single-beam- 
induced backgrounds, the same hadronic event selection 
criteria are applied to separated-beam data, and the num- 
ber of separated-beam events, Ngep, surviving these criteria 
is obtained. The number of beam-associated background 
events, N^g, in the corresponding hadronic event sample 
is given by N^g — fx Ngep, where / is the ratio of the 
products of the pressure at the collision region and the inte- 
grated beam current for colliding and separated-beam runs. 

The beam-associated background can also be determined 
by fitting the distribution of event vertices along the beam 
direction with a Gaussian for real hadronic events and a 
polynomial of degree two for the background, as shown 
in Fig. ^. This was the primary method used for the '99 
data. Figure ^ shows the amount of beam associated back- 
ground versus the scan energy. The differences between R 
values obtained using these two methods to determine the 
beam-associated background range between 0.3 and 2.3%, 
depending on the energy. These differences are included in 
the systematic uncertainty. 

3.3 Hadron Efficiency 

JETSET, the Monte Carlo event generator that is com- 
monly used to simulate e+e^ hadrons, was not intended 
to be applicable to the low energy region, especially that 
below 3 GeV. A special joint effort was made by the Lund 
group and the BES collaboration to develop the LUARLW 
generator, which uses a formalism based on the Lund 
Model Area Law, but without the extreme-high-energy ap- 
proximations used in JETSET's string fragmentation algo- 
rithm [l4|]- The final states simulated in LUARLW 
are exclusive in contrast to JETSET, where they are inclu- 
sive. In addition, LUARLW uses fewer free parameters in 
the fragmentation function than JETSET. Above 3.77 GeV, 
the production of D, D*, Ds, and D* is included in the 
generator according to the Eichten Model JI^ . 

The parameters in LUARLW are tuned to reproduce 14 
distributions of kinematic variables over the entire energy 
region covered by the scan [14, 16|]. For example, the fits 
for parameters tuned at Ecm ~ 2.2 GeV are shown in 
Fig. ||. We find that one set of parameter values is required 
for the CM energy region below open charm threshold, and 
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Figure 3: Distribution of event vertices along the beam line 
(z). A gaussian is used to describe the real beam events, 
while a polynomial of degree two is used to describe beam 
associated background. 
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Figure 4: Amount of beam associated background versus 
scan energy. 



that a second set is required for higher energies. In an alter- 
native approach, the parameter values were tuned point-by- 
point throughout the entire energy range. The detection ef- 
ficiencies determined using individually tuned parameters 
are consistent with those determined with globally tuned 
parameters to within 2%. This difference is included in the 
systematic errors. The detection efficiencies were also de- 
termined using JETSET74 for the energies above 3 GeV. 
The difference between the JETSET74 and LUARLW re- 
sults is about 1%, and is also taken into account in esti- 
mating the systematic uncertainty. Figure ^ (a) shows the 
variation of the detection efficiency as a function of CM 



energy. 
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Figure 5: Comparison of BES data (shaded histogram) and 
LUARLW Monte Carlo data (black histogram) tuned for 
2.2 GeV scan point. 



3.4 Radiative Corrections: (1 + (5) 

The (1 + 6) term is necessary to remove high order effects 
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Different schemes for the radiative corrections were com- 
pared [0 |18|, 0, as reported in Ref. Below 
charm threshold, the four different schemes agree with each 
other to within 1%, while above charm threshold, where 
resonances are important, the agreement is within 1 to 3%. 
The radiative correction used in this analysis is based on 



Ref. [20 1, and the differences with the other schemes are 



included in the systematic error 

3.5 Luminosity Measurement 

The integrated luminosity is determined from 

a e Etrg 

using Bhabha events. The results obtained using 77 and 
dimu events were consistent. The luminosity systematic er- 
ror includes contributions from cut variations, background 
uncertainties, the cross section uncertainty, and the effi- 
ciency uncertainty. The luminosity systematic errors at 
some selected scan energies are shown in Table ^. 

3.6 Trigger Efficiency 

The trigger efficiency was measured by comparing differ- 
ent trigger configurations in special runs at the J/i/'- 
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Figure 6: The cm energy dependence of the detection ef- 
ficiency for hadronic events estimated using the LUARLW 
generator. The error bars are the total systematic errors. 



e^^ = 99.33% 
Chad = 99.76% 
The error on the efficiencies is 0.5 %. 

3. 7 BES Results 

The BES R measurement results | ]l2| , pT| | are shown in Fig. 
/[along with the results from 772, Mark I, and Pluto [ p^ 
23[ p4| ]. Systematic uncertainties are between 6 and 10 % 
and are less than half of the previous uncertainties. The 
average uncertainty is 6.6 %. Tables |l]and^list some of the 
values used in the determination of R and the contributions 
to the uncertainty in the value of i? at a few typical energy 
points in the scanned energy range, respectively. 

4 RECENT CLEO MEASUREMENT 

The CLEO experiment measured R at Ecm = 10.52 GeV, 
just below the T(4S'), and obtained R = 3.56±0.01±0.07 
[p5|]. This very precise measurement (2 %) obtained an 
error similar to that of all previous results in this energy 
region combined, R — 3.579 ± 0.066, corresponding to a 



Table 1 : Some values used in the determination of i? at a few typical energy points. 





jvrobs 
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^had 


+ 


R 


Stat. 
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(%) 






error 


error 


2.000 


1155.4 


19.5 


47.3 


49.50 


1.024 


2.18 


0.07 


0.18 


3.000 


2055.4 


24.3 


135.9 


67.55 


1.038 


2.21 


0.05 


0.11 


4.000 


768.7 


58.0 


48.9 


80.34 


1.055 


3.16 


0.14 


0.15 


4.800 


1215.3 


92.6 


84.4 


86.79 


1.113 


3.66 


0.14 


0.19 



Table 2: Contributions to systematic errors: hadronic selection, luminosity determination, hadronic efficiency determina- 
tion, trigger efficiency, radiative corrections and total systematic error. All errors are in percentages (%). 



EcrniGeV) 


Had. sel. 


L 


Had. eff. 


Trig. 


Rad. corr. 


tot. 


2.000 


7.07 


2.81 


2.62 


0.5 


1.06 


8.13 


3.000 


3.30 


2.30 


2.66 


0.5 


1.32 


5.02 


4.000 


2.64 


2.43 


2.25 


0.5 


1.82 


4.64 


4.800 


3.58 


1.74 


3.05 


0.5 


1.02 


5.14 
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7: A compilation of measurements of R in the cm 
range from 1.4 to 5 GeV. In addition to the BBS 
^ results from 772, Mark I, and Pluto [||, 
are shown. 



1.8 % error [^. Fig. shows the measurements in this 
region prior to CLEO [ 26 1 . 

The CLEO error is smaller than that obtained by BES. 
A breakdown of the error components is given in Table ||, 
along with the corresponding errors at 3.0 GeV from BES. 
Why does CLEO do better? CLEO has much better solid 
angle coverage and a higher detection efficiency, as well 
as a state of the art detector. They also ran with a higher 
luminosity, 1.521 ± 0.015) fb^^, and obtained 4 million 
events. With a much bigger sample, both statistical and 
systematic errors are reduced. 
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Figure 8: R measurements in the 7 to 10 GeV energy re- 



gion. Plot from re f. [26|. 



5 CURRENT STATUS OF R 
MEASUREMENTS 

Burkhardt and Pietrzyk have updated their analysis from 
1995 if] with new results from CMD-2, CLEO, BES, 
and 3rd order QCD for E„n > 12 GeV ||]. They find 
= 128.936 ± 0.046, and A^^^d = 0.02761 ± 
0.00036. In 1995, A^ad = 0.0280 ± 0.0007. The im- 
proved experimental accuracy is primarily due to BES [|[|. 
Fig. ^ shows their current summary of R measurements 
below 10 GeV. The error in the 2 to 5 GeV region is greatly 
reduced because of the BES measurements. Figure |l0| 
shows the breakdown of the error in Aa(Af|) by energy 



(5) 

region. The biggest error contribution to AaJ^g^^ still comes 
from the 1 < Ecm < 5 GeV region! 

Table 3: Comparison of R systematic error contributions 
between BES (at 3.0 GeV) and CLEO (at 10.52 GeV). 
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BES (%) 
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Error (%) 


Error (%) 


Nhad 


3.3 




Backgd./Ev. Modeling 




0.7 
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1.0 


l + S 


1.3 




£had 


2.7 




e-had X (1 + (5) 




1.0 


Sys. 


5.0 


1.8 


Stat. 


2.5 


0.3 


Total 


5.6 


2.0 




Figure 9: Summary of current R values from ref.[|6|]. Com- 
pare with Fig. 0. 



6 FUTURE R MEASUREMENTS 

J. H. Kiihn has used the new BES results to test PQCD, 
but calls for even better precision (2 %) [0]. Certainly, the 
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Figure 10: Contributions to a|^^^^(to^) from the various 
energy regions in magnitude and in uncertainty. Figure 
from ref . [|| . 

precision of R in the 1-2 GeV cm region must be im- 
proved. This is extremely important to improve the preci- 
sion of both and a(M|). The KLOE and PEP-N exper- 
iments are candidates for measuring R in this region. 

As shown in Fig. |l^ the 2-5 GeV region is also very im- 
portant in the determination of Aq;(A/|). The candidates 
for R measurements in this region are BES, PEP-N, and 
CLEO-C. Also improved measurements in this region are 
important to clarify the structure in the charm resonance 
region. The B Factories may also contribute to both energy 
regions using ISR events p7| , |2^ ]. 

The 5 -7 GeV region is also very important. As shown 
in Fig. ^, the R values used by theorists in this region are 
the unpublished Crystal Ball values [^. CLEO-C could 
measure points in this region. 

What is needed to improve R value precision? 

1 . High luminosity large sample 

2. Good soUd angle coverage 

3. Excellent detector with good particle identification 

4. Radiative correction to better than 1 % 



5. More effort on the event generator (LUARLW) 

6. Large sample for tuning the generator 

7. Measure exclusive channels at low energy 

7 ENERGY REACH OF PEP-N 

Although very difficult, extending the energy range to 
above the mass of the ^{2S) should be considered. Cur- 
rently the world's largest sample ( ^4 M events) is that of 
the BES experiment. Some very good physics becomes ac- 
cessible in this region, rr production at threshold can be 
measured to improve the precision of the t mass. Approxi- 
mately one quarter of ■0(25') decays are to Xc and rjc- Run- 
ning on the ip{2S) is a good place to do Xc and rjc physics, 
as well as to study ijj{2S) hadronic decays. The study of 
the hadronic decays might provide a solution to the well 
known pn puzzle Being able to compare J/ip and 

ip{2S) decays is also very beneficial in understanding the 
final states in these decays. Finally it will be very useful 
to have some overlap between PEP-N and CLEO-C as a 
check of systematics in measuring R in this region. 

8 SUMMARY 

Although the precision of R measurements has improved, 
better R precision is still needed to improve the precision of 
the theoretical standard model values of and a{m?^) and 
to test PQCD. There are many interesting possibilities for 
future improvements of R values from KLOE, PEP-N, and 
CLEO-C, as well from the B factories using ISR events. 
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